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I.  INTRODUCTION 

A  numerical  solution  to  the  problem  of  electromagnetic  coupling 
to  a  perfectly  conducting  body  of  revolution  with  an  interior  region 
filled  with  a  loss-free  homogeneous  material  is  developed  here.  Similar, 
but  nevertheless  different,  problems  have  have  already  been  solved.  For 
example,  numerical  solutions  for  homogeneous  material  bodies  of  revolu¬ 
tion  appear  in  [1]  and  [2].  The  problem  of  the  lossy  dielectric  body  of 
revolution  is  addressed  in  [3].  Multilayered  bodies  are  considered  in 
[4]  and  [5].  The  problem  of  electromagnetic  penetration  through  a  sleeve- 
fit  seam  in  an  otherwise  perfectly  conducting  surface  of  revolution  is 
treated  in  [6].  Although  the  seam  in  [6]  may  exhibit  a  dielectric  con¬ 
trast,  the  bulk  of  the  material  inside  the  surface  is  the  same  as  the 
material  outside.  The  concept  of  equivalent  currents  [7,  Sec.  3-5]  is 
used  in  [1 1— [5] .  An  equivalent  aperture  excitation  technique  [8]  and  an 
aperture  loading  technique  are  used  in  [6].  Inspiration  for  the  aperture 
loading  technique  was  drawn  from  [9] . 

Coupling  to  a  perfectly  conducting  body  of  revolution  with  a  homo¬ 
geneous  interior  region  is  considered  in  [10],  but  no  numerical  solution 
is  pursued  there.  Tn  the  present  report,  the  boundary  formulation  con¬ 
sisting  of  [10,  Kqs.  (15)— (18) ]  is  solved  numerically  by  means  of  the 
method  of  moments  [11].  This  formulation  is  only  one  of  several  boundary 
formulations  suggested  in  [10,  Section  2].  More  boundary  formulations 
are  outlined  in  [10,  Sections,  3,  4,  and  5].  Those  in  [10,  Sections  4 
and  5]  take  advantage  of  the  technique  of  subtracting  out  the  short- 
circuit  fields.  This  technique  captures  the  essence  of  the  equivalent 
aperture  excitation  technique.  The  present  report  is  exclusively  de¬ 
voted  to  the  numerical  solution  of  the  particular  boundary  formulation 
consisting  of  [10,  Eqs.  (15)-(18)]. 

II.  I'OKMUl.ATION  OF  Till'  I’ROliLEM 

The  problem  is  that  of  electromagnetic  coupling  to  the  perfectly 
conducting  body  of  revolution  with  homogeneous  interior  region  shown  in 
Fig.  1.  Figure  1  shows  a  cross  section  of  the  configuration  in  a  plane 
containing  the  axis  (z  axis)  of  the  conducting  body  of  revolution.  This 
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Fig.  1.  Original  problem. 
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conducting  body  is  crosshatched .  Homogeneous  material  with  permeability 
y+  and  permittivity  £+  fills  the  exterior  region  V+  bounded  by  the  surface 
S+  and  the  aperture  A.  What  is  called  S+  is  not  a  single  surface,  but  two 
surfaces,  one  below  A  and  one  above  A.  The  interior  region  V  bounded  by 
the  surface  S  and  the  aperture  A  is  filled  with  homogeneous  material 
characterized  by  y  and  e  .  All  three  surfaces  S+,  S  ,  and  A  are  surfaces 
of  revolution.  The  electromagnetic  excitation  consists  of  impressed  electric 
and  magnetic  current  sources  and  M^+  in  V+  and  sources  J*  and  in  V  . 

+  4*  + 

The  resulting  electric  and  magnetic  fields  in  V  are  called  E  and  H  .  The 
fields  in  V  are  called  _E  and  H  .  The  fields  E~  and  H’  are  unknown. 
Obviously,  solution  for  them  constitutes  a  complete  solution  to  the  problem 
presented  in  Fig.  1.  However,  our  solution  will  be  for  equivalent  surface 

+  4  4 

currents  on  S  and  A  rather  than  for  the  fields  E  and  H  .  It  is  possible 

+■ 

to  construct  E  and  H~  from  these  equivalent  surface  currents,  but  this  is 
not  done  in  the  present  report. 

The  equivalence  principle  [7,  Sec.  3-5]  is  used  to  obtain  the  situ- 

i+  i+ 

ations  shown  in  Figs.  2  and  3.  In  Fig.  2,  (J_  ,  M  )  and  equivalent  cur¬ 

rents  J+,  .J,  and  M  radiate  in  the  presence  of  the  homogeneous  medium 
(y+,  e+)  to  produce  (E+,  H+ )  in  V+  and  zero  field  elsewhere.  Here,  J+  is 
an  electric  current  on  S+,  J  is  an  electric  current  on  A,  and  M  is  a  mag¬ 
netic  current  on  A.  In  Fig.  3,  (J1  ,  M1  )  and  equivalent  currents  -J  ,  -J, 
and  -M  radiate  in  the  presence  of  the  homogeneous  medium  (y  ,  £  )  to  pro¬ 
duce  (E  ,  H  )  in  V  and  zero  field  elsewhere.  Here,  -J  is  an  electric  cur¬ 
rent  on  S  ,  -J  is  an  electric  current  on  A,  and  -M  is  a  magnetic  current 
on  A.  Since  the  surfaces  S+  and  S  are  perfectly  conducting  in  the  original 
problem  of  Fig.  1,  only  equivalent  electric  currents  are  needed  on  them  in 
Figs.  2  and  3.  The  choice  of  -J  rather  than  J  on  S  in  Fig.  3  is  due  to 
personal  preference.  However,  the  minus  sign  relationship  between  the 
aperture  currents  in  Figs.  2  and  3  is  mandated  by  the  zero  field  stipula¬ 
tions  in  Figs.  2  and  3  and  continuity  of  the  tangential  components  of  the 
electric  and  magnetic  fields  across  the  aperture  in  Fig.  1.  If  the  zero 
field  stipulations  in  Figs.  2  and  3  are  enforced,  then  the  minus  sign 
relationship  between  the  aperture  currents  in  Figs.  2  and  3  ensures  con¬ 
tinuity  of  tangential  fields  across  the  aperture  in  Fig.  1.  However, 
if  these  zero  field  stipulations  are  not  enforced,  then  the  minus  sign 
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relationship  between  the  aperture  currents  does  not  ensure  continuity 
of  the  tangential  fields  across  the  aperture  in  Fig.  1. 

The  equivalence  principle  states  that  the  equivalent  currents  in 
Figs.  2  and  3  are  unique,  but  does  not  tell  how  to  determine  them.  The 
expressions  in  [7,  Sec.  3-5]  for  the  equivalent  currents  in  terms  of  the 
tangential  components  of  the  fields  can  not  be  used  because  the  fields 
are  not  known  either.  We  determine  the  equivalent  currents  by  enforcing 
the  boundary  conditions  for  the  fields  in  Fig.  1.  These  boundary  con¬ 
ditions  are: 

1)  The  tangential  components  of  the  electric  field  vanish 
on  the  conducting  surface  S+. 

2)  The  tangential  components  of  the  electric  field  vanish  on 
the  conducting  surface  S  . 

3)  The  tangential  components  of  the  electric  field  are  con¬ 
tinuous  across  the  aperture  A. 

4)  The  tangential  components  of  the  magnetic  field  are  con¬ 
tinuous  across  the  aperture  A. 

The  equations  for  these  four  boundary  conditions  are 


—  E+  =0 
r  +  -tan 


—  E~ 

+  —tan 

n 


J_  E+ 

+  —tan 


n 


i_  E- 

+  —tan 


!L  =  h 

—tan  —tan 


on  S 


on  S 


on  A 


on  A 


(la) 

(lb) 

(lc) 

(ld) 


The  subscript  "tan"  denotes  tangential  components  on  the  surface  in  ques¬ 
tion.  The  factor 


v/4 


JL_ 

n  u 

was  included  in  (1)  to  make  (la),  (lb),  and  (lc)  compatible  with  (Id), 


(2) 
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In  view  of  Figs.  2  and  3,  (1)  can  be  rewritten  as 


E,  (J  +J,  M)  =  ~  E~ 

+  tan  —  —  —  +  tan 


+  «■  ">  ‘ 


-  i  E1* 
+  —tan 

n 


-  -t:  £.n«.+  +i-  "> 


~  E  (J  +J,  M)  =  ET 
+  —tan  —  —  —  +  —tan 

n  n 


.  JL  Ei+  -  -L  Ei- 

+  —tan  +  —tan 


H+  (J+  +J,  M)  -  H~  (J~  +J,  M)  =  H*+  -  H*' 
—tan  —  —  —  — tan  —  —  —  —tan  —tan 


on  A  (3d) 


where  E  (J,  M)  denotes  the  tangential  components  of  the  electric  field 
— tan  —  — 

due  to  the  electric  current  J  and  the  magnetic  current  M,  both  radiating 

+  +  + 

in  the  medium  characterized  by  (u  ,  e~) .  Similarly,  HT  (J,  M)  denotes 

-tan  -  ~  i+ 

the  tangential  components  of  the  corresponding  magnetic  field.  E  and 

i±  i±~Cani± 

H  are  the  tangential  electric  and  magnetic  fields  due  to  (J  ,  M  ) 
tan  +  +  jf. 

radiating  in  the  medium  characterized  by  (y“,  e“) .  Since  J  and  are 
disjoint,  the  +  sign  between  them  in  the  expression  (J*  +J)  means  union 
rather  than  simple  addition.  Similarly,  (J  +J)  denotes  the  union  of  J 
with  J.  It  does  not  matter  whether  the  tangential  electric  fields  on  the 

left-hand  sides  of  (3a)  and  (3b)  are  evaluated  just  inside  S*  or  just 

+  +  f 
outside  S  because  they  are  continuous  across  S  .  The  Ampere  s  law  con¬ 
tributions  cancel  out  of  the  left-hand  sides  of  (3c)  and  (3d).  By  Ampere's 
law  contribution  we  mean  the  contribution  to  the  field  due  to  the  value  of 
the  current  at  the  field  point.  The  qualification  "on  A"  in  (3c)  and  (3d) 
means  not  only  evaluation  on  A  but  also  suppression  of  all  Ampere's  law 
contributions.  The  set  of  equations  (3)  is  essentially  the  same  as  [10, 

Eqs.  (15) — (18) ] .  It  is  shown  in  [10]  that  these  equations  suffice  to 

+ 

uniquely  determine  the  equivalent  currents  J  ,  and  M. 

The  four  field  operators  which  give  the  electric  and  magnetic  fields 
due  to  electric  and  magnetic  currents  appear  in  (3).  However,  it  is  evident 
from  [7,  Eq.  3-79)]  that 


(0,  M>  -  -  nV  (i  K.  0) 


HZ  (°»  M) 
—tan 


7~2  i£an<7  0) 

(rf)  n 


(5) 


where 


The  seemingly  superfluous  pair  of  factors  H+  and  l/n+  on  the  right-hand 
sides  of  (4)  and  (5)  were  included  in  order  to  liken  the  first  arguments 

*  t 

of  E  and  H  to  electric  currents.  In  view  of  linearity  of  the  field 
—tan  —tan  J 

operators,  substitution  of  (4)  and  (5)  into  (3)  gives 


~  F+  (J++J,0)  +  H+  (~~  M,  0)  =  —  E1+ 
+  —tan  -  —  -tan  +  —  +  —tan 

n  n  n 


<t:  h.o> 


~  Ei_ 

H - tan 


-  -+  <a„y>j.o)  -  J-  ♦  a;an(i  «,»>  +  «;an(i »,<» 

n  n  n  n 

=  ■— -  E l+  -  ~  E1”  on  A  (7c) 

n+  -tan  n+  -tan 

-  H*  (J++J,0)  -  H~  (j“+j,0)  -  ~  E+  (~  M,0) - E~  ( ~  M,0) 

-tan  -  -  -tan  -  —  +  —tan  '  +  -  ,  -  .  2  -tan  +  - 

n  n  (n  )  n 

=  ui+  -  H i_ 

— tan  —tan  on  A  (7d) 

Unlike  the  left-hand  sides  of  (3)  which  contain  all  four  field  operators, 
the  left-hand  sides  of  (7)  contain  only  the  operators  which  give  the  electric 
and  magnetic  fields  due  to  electric  currents.  Equations  (7)  are  solved  by 
means  of  the  method  of  moments  in  section  III. 


III.  METHOD  OF  MOMENTS  SOLUTION 


The  method  of  moments  solution  to  (7)  is  obtained  bv  first  express¬ 
ing  the  unknown  currents  (J++J),  (.J  +  J_) ,  and  M  as  linear  combinations  of 
known  expansion  functions,  then  taking  the  symmetric  product  of  (7)  with 
each  member  of  a  set  of  testing  functions,  and  finally  solving  the  resulting 
set  of  simultaneous  equations. 
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Momentarily  ignoring  the  fact  that  (J_  +J)  and  (J  +J)  are  supposed 
to  be  equal  on  A,  we  write 


N+-2 

l 

It+ 

t+  N+_1 

J  +  1 

on  (S++A) 

(8) 

J=i 

nj 

-  nj 

nj  -nj 

N~-2 

N~-2 

i 

I 

1*7 

+  l 

i  . 
% 
t 

on  (S  +A) 

(9) 

j  =  l 

nj 

-j  j=i 

nj  -nj 

J++J  =  }' 


l+l  =  l 

n 


where  (S  +A)  is  the  surface  which  is  the  union  of  S  and  A.  Similarly, 
(S  +A)  is  the  union  of  S  and  A.  Since  S+,  S  ,  and  A  are  surfaces  of 
revolution,  (S++A)  and  (S  +A)  are  also  surfaces  of  revolution.  In  (8) 

t  +  (f>  f  £  + 

and  (9),  I  .  and  I  are  unknown  coefficients.  The  vectors  J  .  and 


nj 


nj 


-nj 


J^j  in  (8)  and  (9)  are  expansion  functions  defined  by 


jt;  =  u 
a  j  -t 


t;  (f)  .  . 

-j _ eJn<J> 

P 


.  !ilil 


|j  -  1,2....  IC-2 

In  -  0,  ‘1,  *2,... 

|j  =  1,2,  . . .  (N+-l  or  N~-2) 

) n  =  0,  +1,  f 2, . . . 


(10a) 


(10b) 


The  right-hand  sides  of  (10)  were  obtained  by  attaching  the  superscript  ± 
to  various  quantities  in  the  version  of  [12,  Eqs.  (2)  and  (3)]  with  P 

replaced  by  X.  The  notation  ±  in  (10)  denotes  quantities  defined  on 

;  ±  + 

(S  +A) .  The  subscript  j  which  runs  from  1  to  either  N  -2,  N  -1,  or 

N  -2  in  (10)  is  not  to  be  confused  with  the  j  which  appears  in  the  argu¬ 
ment  of  the  exponential  in  (10).  The  latter  j  is  /-l. 


The  quantities  decorated  with  the  notation  ±  on  the  right-hand 

■f 

side  of  (10)  depend  on  the  generating  curve  of  the  surface  (S_+A) .  The 
generating  curve  of  a  surface  of  revolution  is  the  plane  curve,  which  when 
rotated  about  the  z  axis,  generates  the  surface  of  revolution.  For  com- 

4* 

putations,  the  generating  curve  of  (S~+A)  is  approximated  by  choosing  a 

-■*  +  + 
succession  of  points  t.,  j=l,2,...N~,  on  the  generating  curve  of  (S  +A) 

3  _ 

and  then  connecting  them  with  straight  line  segments.  The  point  tj  is 

M  + 

called  a  data  point.  Its  location  is  specified  by  its  distance  p~  from 

—  1  J 

the  z  axis  and  its  z  coordinate  z..  Hence,  we  write 
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Location  of  tj  =  (fy  Zj),  j  =  1,2,... (11) 

4.  —  +  —  +■  —  + 

The  generating  curve  of  (S  +A)  starts  at  t",  goes  to  t-,  then  to  t . , 

—  ±  ±  £  J 

and  so  forth  until  the  last  data  point  t  +  is  reached.  A  typical  arrange- 

N"  + 

ment  of  data  points  on  the  generating  curve  of  (S  +A)  is  shown  in  Fig.  U. 
Figure  5  shows  an  arrangement  of  data  points  on  the  generating  curve  of 
(S  +A) .  The  first  and  last  data  points  on  the  generating  curve  of  (S++A) 
are  on  the  z  axis.  The  last  two  data  points  on  the  generating  curve  of 
(S  +A)  overlap  with  the  first  two.  All  data  points  on  A  are  common  to 

■j*  — 

the  generating  curves  of  both  (S  +A)  and  (S  +A) .  There  are  M  data  points 
on  A.  The  first  data  point  on  A  is  the  point  where  A  meets  the  lower 
part  of  S  .  This  data  point  is  the  (M  +l)th  data  point  on  the  generating 

4- 

curve  of  (S  +A).  The  last  data  point  on  A  is  the  point  where  A  meets  the 
upper  part  of  S+. 


In  (10),  0  is  the  angular  distance  from  the  positive  x  axis. 

<t>  is  measured  in  the  xy  plane  and  toward  the  positive  y  axis.  In  (10b), 
u.  is  the  unit  vector  in  the  <J)  direction.  No  ±  needs  to  be  attached  to 

-v  + 

<)>  and  u^  because  they  do  not  depend  on  the  generating  curve  of  (S  +A)  . 

+ 

Henceforth,  the  expression  "generating  curve  of  (S  +A)"  will 
mean  the  series  of  straight  line  segments  connecting  the  data  points 

_4  4  4 

tV,  j=l,2,...N‘.  In  (10),  t‘  is  the  arc  length  along  the  generating 
1  +  + 
curve  of  (S  +A)  and  u  +  is  the  unit  vector  in  the  t  direction.  We 

+  ~ ^  _+ 

choose  t  to  be  zero  at  the  data  point  t.  where  the  generating  curve  of 
1  _+ 

(S  +A)  starts.  Previously,  only  the  location  of  t.  was  specified.  The 
-  *  +  3  _+ 

value  of  t.  is  now  defined  to  be  the  value  of  t  at  the  location  of  t  . 

1  +  f  4  .4  1 

In  (10a),  T~  (t-)  is  the  triangle  function  which  starts  at  t"  =  y, 

’  +_4  4J_4 

reaches  the  maximum  height  of  unity  at  t  =  t,  . ,  and  ends  at  t_  =  t.  ». 

4  4  j  +  ^  4  _4  j  +  ^ 

In  (10b),  P,(t  )  is  the  unit  pulse  function  which  starts  at  t  =  t. 

j  4  _4  4  4  J 

and  ends  at  t~  =  tj+^*  Naturally,  the  functions  Tj  and  Pj  are  real.  They 

are  essentially  the  same  as  the  functions  T,  and  P,  shown  in  [12,  Figs.  1 
-  ,  3-4  3 

and  2].  The  tj  used  in  [12]  corresponds  to  tj  of  the  present  report.  The 

notation  of  the  present  report  was  obtained  by  replacing  the  superscript  - 

on  t.  in  [12]  by  a  horizontal  bar  directly  above  t,  and  then  appending 
3  *  3  4 

the  superscript  f.  In  (10a),  p  is  the  distance  from  the  z  axis  at  t  . 
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*  .  4  4 

In  (10b),  p.  is  the  value  of  p  at  the  center  of  the  pulse  function  P., 


Because  of  the  nature  of  the  approximate  generating  curve  of  (S  +A) , 

> 


’  i  <pj  * 


(12) 


-4  _  + 

where,  as  previously  defined,  p.  is  the  distance  of  the  data  point  t. 

4  4*4  J 

from  the  z  axis.  With  the  advent  of  t~  and  p ",  p7  becomes  the  value 

4  4-4  d 

of  p  at  t "  =  t ' . 

As  is  evident  from  the  range  of  values  of  j  in  (10a),  a  peak  of 
a  triangle  function  is  placed  at  each  data  point  on  the  generating  curve 

— - j-  — 

of  (S  +A)  except  the  first  data  point  t,  and  the  last  data  point  t 

-+  1  _f  N+ 

As  shown  in  Fig.  A,  the  data  points  t,  and  t  .  lie  on  the  z  axis.  No 

1  _+ 

peak  of  a  triangle  function  is  needed  at  t^  because  the  function 
T\(t+)/p+  in  (10a),  being  equal  to  1/pt  for  tl"  £  t+  £  t^,  is  just  as 
large  at  t.  as  it  is  at  t..  Similarly,  no  peak  of  a  triangle  function 

— H  ^  +  +  + 

is  needed  at  t  ,  because  the  function  T  (t  )/p  in  (10a)  holds  its 
N  +■  -+  N+-2  +  _+ 

value  I/p  at  t  =  t  .  all  the  way  out  to  t  =  t  . 

N+-1  N  -1  N+ 


Tt  was  necessary  to  overlap  the  last  two  data  points  with  the 
first  two  in  Fig.  5  in  order  to  induce  a  computational  procedure  which 
places  no  peak  of  a  triangle  function  at  the  first  and  last  data  points 
to  obtain  a  peak  of  a  triangle  function  at  every  data  point  on  the 
generating  curve  of  (S  +A) .  The  above  mentioned  computational  procedure 
was  originally  designed  for  a  generating  curve  such  as  that  of  (S++A) 
whose  first  and  last  data  points  are  on  the  z  axis  where  no  peaks  of 
triangle  functions  are  wanted.  The  second  summation  with  respect  to  j 

in  (9)  is  terminated  at  j  =  N  -2  rather  than  at  j  =  N  -1  because,"  due  to 

i- 

the  overlapping  data  points  in  Fig.  5,  the  vector  function  J  ,  de- 

A_  “n’N  _1 

fined  by  setting  j  =  N  -1  in  (10b)  is  a  duplicate  of  J  . 

~ nl 

Because  each  triangle  function  in  (10a)  is  continuous  in  t  , 

♦  + 
the  t  components  of  (8)  and  (9)  are  also  continuous  in  t  .  As  a  result, 

(8)  and  (9)  do  not  give  rise  to  any  circular  loops  of  filamentary  electric 

charge.  Any  circular  loop  of  such  charge  is  not  physically  realizable 

because  an  infinite  amount  of  energy  is  required  to  assemble  it.  Since 

each  pulse  function  in  (10b)  begins  and  ends  abruptly,  d iscont inui t ies 
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in  the  <p  components  of  (8)  and  (9)  are  allowed  at  all  the  data  points. 

+ 

Unlike  discontinuities  in  the  t  component  of  electric  current,  dis¬ 
continuities  in  the  <}>  component  of  electric  current  are  not  accompanied 
by  concentrations  of  electric  change. 

Since  the  data  points  of  the  generating  curves  of  (S++A)  and  (S  +A) 

coincide  with  each  other  on  A,  the  vector  functions  Jt.  and  J  .  on  A  in 

— nj  -nj 

(8)  are  the  same  as  those  on  A  in  (9).  Hence,  (8)  and  (9)  can  be  made 
equal  on  A  by  equating  the  coefficients  of  these  vector  functions.  We 
write 


j++j  =  l  j  it+jt+  +  l  i^+.  +  L  +  L 

*  "  n  j-1  nj_nj  j-1  nJ-nJ  j  =M  +M  nj^j  j  =M  +M  nj“ 


<t>+  t<k 

nj-nj 


+  y  i  i  ,j  . ,  .  +  i  i  .  .r 

L  ,L  nj— n,  j+fr-l  n,  j+M-n,  j+M+ 

n  j=l  j-i 


1  +1  =  n  j-1  Iflj ^  +  j-1  +  j  =M~+M  +  j  Jr+M1"™ 


M  M-l  , 

r->  [  “  r*  © — 

+  £  Inj^n,j+M  -1  +  In,j+M 


,t+  Tt-  . 

4m+  =  4m-  on  A 


(15a) 


Jt+  -  jt_ 

^n.j+K^-l  — n, j+M“-l ’ 


=  2,3, ...M-l 


(15b) 


(15c) 


•  j-1-2--"-1  <i5« 

The  coefficients  that  were  equated  in  (8)  and  (9)  have  been  relabeled 

I  .  and  I  ....  in  (13)  and  (14).  The  first  )  in  (13)  is  a  vector  func- 
nj  n, j+M  ft 

tion  confined  to  S+.  The  vector  functions  J^+  anc*  in  0-3) 

straddle  S+  and  A.  Except  for  the  overflow  of  J^f  and  ^ 


onto 


13 


S+,  the  second  £  In  (13)  is  a  vector  function  confined  to  A.  The  first  £ 


in  (14)  is  confined  to  S  .  The  vector  functions  J*  and  . .  , 

~nM  — n,MiM  — l 


n 

in  (14) 


straddle  S  and  A.  Except  for  the  overflow  of  and  JC  . . , w-  ,  - 

— nM  — n ,M+M  -1  onto  S  , 

the  second  £  in  (14)  is  confined  to  A.  The  right-hand  sides  of  (13)  and 
(14)  are  equal  on  A  because  (15)  implies  that  the  portion  on  A  of  the 
second  £  in  (13)  is  equal  to  the  portion  on  A  of  the  second  £  in  (14). 


To  condense  the  notation  in  (13)  and  (14),  we  define  expansion 


1  +  2  + 

functions  J  .  and  J  .  by 

-nj  — nj 


j  =  1,2,.. .M~-l 
j  =  M±,  M±+l, . . .  2M±-1 
j  =  2M+ ,  2M*+1, . . ■N±+M±-M-2 


^n, j-N1+2M+l 


,  j  =  N~+M+-M-l . . . 


2N+-2M-2 


or  2M  -2M-3 


(16) 


j2: 

-nj 


and  coefficients  I  .  by 
nj  J- 


^n,j+M±-l 


,j+M±-M 


Xnj 


nj 


j  =  1,2, ...M 


j  =  M+l,  M+2, . . .2M-1 


,  j  =  l,2,..tT-l 


n,j-M+-+l 

- 

n.j-M^M 
n,j-Ni+2M+l’ 


1  =  M" ,  M+l ,  . . .2M  -1 


,  j  =  2M“ ,  2M++l,...N±+M±-M-2 


j  =  fT+M*-M-l, . • . 


\2N+-2M-2 


(17) 


(18) 


|or  2N  -2M-3 

The  new  vector  functions  (16)  and  (17)  and  the  new  coefficients  (18)  re¬ 
duce  (13)  and  (14)  to 
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r  + 

2N  -2M-2 

1-1-  1-1- 

2M-1 

J  + J  = 

l 

l 

I1+J1+ 

+ 

l 

i  .J2+ 

(19) 

n 

j=1 

nj-nj 

j-1 

nj-ni 

2N_-2M-3 

l 

2  M-l 

J-+J  = 

l 

t1Tj1T 

+ 

V 

1  ^j2~ 

(20) 

n 

j-1 

nj-nj 

j-l 

nj-nj 

The  first  J  in  (19)  is  on  S+.  The  functions  and  .  in  (19) 

j  — nl  — n,2M-l 

+  2+  2+ 
straddle  S  and  A.  Except  for  the  overflow  of  J  ,  and  J  ...  ,  onto 
+  — nl  — n,2M-l 

S  ,  the  second  £  in  (19)  is  confined  to  A.  The  first  £  in  (20)  is 

j  j 

2-2- 

on  S  .  The  functions  J  ,  and  J  ...  ,  in  (20)  straddle  S  and  A.  Except 

— nl  — n, 2M-1  v  r 

2“  2—  — 

for  the  overflow  of  J  ,  and  J  ,  onto  S  ,  the  second  T  in  (20)  is  con- 
— nl  —n,2M-l  t 

fined  to  A.  The  right-hand  sides  of  (19)  and  (20)  are  equal  to  each  other 
on  A. 


The  equivalent  magnetic  current  M  in  (7)  is  expanded  as 


.  2M-3 

M  =  n  I  l  v  J  "T  (21) 

n  j-1  nJ^ 

3+ 

where  V  .  is  an  unknown  coefficient  and  J  .  is  an  expansion  function 
nj  ~nj 

given  by 


.  j  =  1,2, ...M-2 


j<H 

^n,j+ttMl+2  ,  j  =  M-l,  M....2M-3 


(22) 


The  t+  component  of  (21)  is  continuous  so  that  no  circular  loops  of  fila¬ 
mentary  magnetic  charge  appear  in  the  aperture  in  Figs.  2  and  3.  The  4> 
component  of  (21)  can  have  a  discontinuity  at  each  data  point  in  the  aper¬ 
ture. 


Now  that  the  expansions  ( 19)  —  ( 21 )  have  been  written  for  the  equi¬ 
valent  currents,  the  next  step  in  the  method  of  moments  is  to  take  the 
symmetric  product  of  (7)  with  each  member  of  a  set  of  testing  functions. 
The  symmetric  product  of  two  vectors  W  and  E  is  called  <W,E>  and  is 
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defined  to  be  the  integral  of  the  dot  product  of  W  with  JE  over  S  ,  S  , 
and  A.  First,  (7a)  is  tested  with  the  vector  function  W^T  defined  by 

1+  1+  *  + 

W  .  =  (J  ,)  ,  i  =  1,2, ...2N  -2M-2  (; 

ni  — ni 

where  the  asterisk  denotes  complex  conjugate.  Then,  (7b)  is  tested 


W1"  -  (J1!)*  ,  i  =  1,2, . . .2N  -2M-3 
— ni  — ni 


For  the  sake  of  argument,  we  write 


E*  (J++J,0)  +  H*  (~r  M,0)  =  \  Ef+ 

+  -tan  —  —  —tan  +  —  +  —tan 

n  n  n 


on  (S  +A) 


T:  s;,n(i'w.0)  +  h.o)  -  -  ±  on  (s'+a) 


2+  2- 

Testing  functions  and  are  defined  by 


w2;  =  (J2+)* 

— ni  — ni 


i  =  1,2,... 2M-1 


2-  2-  * 

-  (4i> 


i  =  1,2,... 2M-1 


For  a  particular  value  of  i,  (25)  is  tested  with  W  , ,  (26)  is  tested 

2-  _nl 
with  W^,  and  the  results  are  added  to  obtain  one  equation.  Letting  i 

run  from  1  to  2M-1,  we  obtain  2M-1  equations.  This  kind  of  testing  is 

equivalent  to 

2+  + 

1)  Testing  (7a)  with  the  part  of  W  .  on  S  , 

— ni 

2- 

testing  (7b)  with  the  part  of  on  S  , 

2+ 

testing  (7c)  with  the  part  of  W  ,  on  A, 

— ni 

and  adding  the  results  to  obtain  one  equation. 


2)  Testing  (7c)  with  W  7  for  i  =  2,3,...M-1 

— ni 


2+  + 

3)  Testing  (7a)  with  the  part  of  W  on  S  , 

— nM 

2- 

testing  (7b)  with  the  part  of  W  on  S  , 

nM 

24- 

test  ing  (7c)  with  the  part  of  W  on  A, 

“nM 

and  adding  the  results  to  obtain  one  more  equation. 

4)  Testing  (7c)  with  W2"^  for  i  =  M+l,  M4-2  , .  .  .  2M-1 . 

— ni 

Statement  1)  is  justified  in  the  following  manner.  Testing  (25)  with 
2+  + 

the  part  of  on  S  is  equivalent  to  testing  (7a)  with  the  part  of 

24-  +  + 

W  ,  on  S  because  (7a)  and  (25)  are  identical  on  S  .  Likewise,  testing 

_nl  2-  - 

(26)  with  the  part  of  W  on  S  is  equivalent  to  testing  (7b)  with  the 
2  ^  1 

part  of  W  on  S  .  It  is  evident  from  (27),  (28),  (17),  and  (15)  that 
— nl 

w2+  =  w27  on  A,  i  =  1,2,... 2M-1  (29a) 

— ni  — ni 

24- 

Hence,  testing  (25)  with  the  part  of  on  A,  testing  (26)  with  the 

part  of  W2  on  A,  and  adding  the  results  is  equivalent  to  testing 

_Tli  24- 

the  sum  of  (25)  and  (26)  with  the  part  of  W  .  on  A.  Because  the  sum 

— nl 

of  (25)  and  (26)  on  A  is  precisely  (7c),  this  amounts  to  testing  (7c) 

24- 

with  the  part  of  on  A.  Consequently,  statement  1)  is  true. 

Since  all  the  testing  functions  in  statements  2)  and  4)  are 
confined  to  A,  (29a)  implies  that 


24-  2- 

W  .  =  W  . 
— ni  — ni 


for  all  the  values  of  i  in  statements  2)  and  4).  Hence,  for  these 

24-  2- 

values  of  i,  testing  (25)  with  W  testing  (26)  with  W  .,  and  adding 

— nl  — nl 

24- 

the  results  amounts  to  testing  (7c)  with  because  the  sum  of  (25) 
and  (26)  on  A  is  (7c).  Thus,  statements  2)  and  4)  are  true.  Justifi¬ 
cation  for  statement  3)  is  similar  to  that  for  statement  1)  and  is 
left  as  an  exercise  for  the  interested  reader. 


Finally,  (7d)  is  tested  with 
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Substituting  the  expansions  ( 1 9>  —  (21 )  into  (7)  and  then  test¬ 
ing  (7)  with  (23),  (24),  (27),  (28),  and  (30)  in  the  previously 
described  manner,  we  obtain  the  matrix  equations 


T  X  =  B  ,  n  =  0,  +1,  ±2,  . .  . 
n  n  n 


(31) 


Here,  X^  is  a  column  vector  of  the  unknown  coefficients  in  (19)-(21). 

Also,  T  is  a  square  matrix  called  the  moment  matrix.  Finally,  B  is 
n  n 

a  column  vector  callc=a  the  excitation  vector.  The  elements  of  B  are 

n 

obtained  by  testing  .  e  known  fields  Ejf1  and  on  the  right-hand 


—tan 


—tan 


sides  of  (7) . 


The  column  ve.,  tor  X  is  defined  by 
n  J 


X  =  [i1+ 

n  n 


-1- 


V 


(32) 


where  the  tilde  (~)  denotes  the  transpose  of  a  column  vector  and  where 
the  jth  elc 
I  .,  and  V 


the  jth  elements  of  I3+,  I3  ,  I  ,  and  V  are,  respectively,  I^T,  I3., 

n  n  n  n  nj  nj 


nj’ 


*  .  • 
nj 


The  moment  matrix  T^  is  given  by 


z11+ 

n 

0 

z12+ 

n 

-y13+ 

n 

0 

11- 
n  z 

r  n 

\<2~ 

-Y13- 

n 

z21+ 

n 

21- 

nrZn 

z22+  +  n  z22- 

n  r  n 

_y23+  _  y23- 
n  n 

y31+ 

n 

Y31- 

n 

32+  32- 

Y  +  Y 

n  n 

z33++(— )z33- 
n  n  n 

_ r _ — 

(33) 


where 


n 


(34) 


The  ijth  elements  of  the  superscripted  Z  and  Y  submatrices  in  (33) 

n  n 

are  given  by 


1 
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zpq±  =  _<wp±  1  0)> 

nij  ni  *  ±  E  ’  0) 


P  =  1,2,3 


YPqN-<WP:,  H+(JqT  ,  0)> 
nij  ni  — nj 


q  =  1,2,3 


3-  3- 

The  new  vector  functions  J  .  and  W  .  unavoidably  introduced  in 

— nj  — ni 

(35)  and  (36)  are  defined  by 


(35) 

(36) 


(37) 


4-  + 

No  superscript  "tan"  is  needed  on  the  field  operators  E~  and  H_  in  (35) 

pi 

and  (36)  because  the  testing  function  W  .  is  tangential  to  the  surface 
.  I  — ni 

(S  +A).  Because  W'  and  J  .  have  the  dimension  of  (1/distance)  and 
— ni  — nj 

1  +  i 

because  the  operators  — r  E  and  H  are  dimensionless,  expressions  (35) 
and  (36)  are  dimensionless.  Hence, the  elements  of  the  moment  matrix  T 

n 

are  also  dimensionless. 


-4 

The  excitation  vector  in  (31)  is  given  by 


B  =  [Vll+  Vil_  V1  I1]  (38) 

n  n  n  n  n 


where  the  ith  elements  of  V^+,  V*3  ,  V*  and  I*  are  given  by 

n  n  n  n 


V 

V 

V 

I 


il+ 

ni 

il- 

ni 

i 

ni 

i 

ni 


<W1+  ,  i  E1+> 
nx  + 

n 

-  <»;i .  ^ 

^  -  <£  •  ^ 
n  n 


<W3*  ,  Hi+>  -  <W3*  H1-> 

ni  ni 


Ei-> 


(39a) 

(39b) 

(39c) 

(39d) 


The  superscript  i  in  (38)  and  (39)  is  not  to  be  confused  with  the  sub- 

i  + 

script  i  in  (39).  The  former  i  indicates  dependence  on  the  fields  E 
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and  H1  .  The  latter  i  denotes  the  ith  element  of  the  column  vector 

i  -  i  — 

in  question.  No  superscript  "tan"  is  needed  on  _E  and  H  in  (39) 
because  all  the  testing  functions  in  (39)  are  tangential  vectors. 

In  general,  the  method  of  moments  gives  one  matrix  equation  in 
which  all  the  unknowns  are  coupled  to  each  other.  However,  in  the 

— ►  -V 

matrix  equations  (31)  there  is  no  coupling  between  X  and  X  for  m  ^  n. 

n  m 

This  lack  of  coupling  is  a  consequence  of  the  nature  of  the  field 
operators  in  (7),  the  e''n’5  dependence  of  the  expansion  functions,  the 

—  i  ft'* 

e  dependence  of  the  testing  functions,  and  the  fact  that  the 

integral  from  0  to  2~  of  the  product  of  e  with  e^n<^  is  zero  whenever 

m  ^  n.  The  two  field  operators  which  give  the  tangential  electric 

and  tangential  magnetic  fields  due  to  an  electric  current  appear  in  (7). 

An  e)nv  electric  current  produces  an  e^n<^  electric  field  [13,  Appendix  B] 
i  nv 

and  an  eJ  v  magnetic  field  [13,  Eq.  (9)].  In  the  previous  sentence,  the 

i n6  i 

modifier  eJ  means  that  the  tp  dependence  is  described  by  eJ  v.  Actually, 

[13,  Eq.  (9)[  represents  the  operator  ~n  *  H  rather  than  ~H  where  H  de¬ 
notes  the  magnetic  field  due  to  an  electric  current  and  n  is  the  unit 
normal  vector  which  points  outward  from  the  surface  of  the  body  of  revo¬ 
lution  in  [13,  Fig.  1].  However,  the  effect  of  n  in  n  *  H  is  to  discard 
the  normal  component  of  H  and  to  rotate  the  tangential  component  of  H 
through  an  angle  of  90°  so  that  if  n  *  H  is  ejn<^,  then  the  tangential 
component  of  II  is  also  e''n<^. 

Equation  (7d)  can  be  called  a  magnetic  field  equation  because  it 
comes  from  the  boundary  condition  that  the  tangential  components  of  the 
magnetic  field  are  continuous  across  the  aperture.  Since  (7d)  is  a 
magnetic  field  equation,  we  were  tempted  to  multiply  it  by  n x  in  order 
to  make  it  conform  to  [13,  Eq.  (3)]  and  [14,  Eq.  (1)].  The  nx  is  vital 
to  [14,  Eq.  (1)].  If  the  nx  were  removed  from  [14,  Eq.  (1)],  its  method 
of  moments  solution  using  the  testing  and  expansion  functions  defined  in 
[14]  would  fail  for  any  conducting  body  of  revolution  with  a  symmetric 
generating  curve.  The  conducting  sphere  in  [14,  Fig.  1]  is  an  example 
of  a  conducting  body  of  revolution  with  a  symmetric  generating  curve. 

The  above  mentioned  failure  of  the  method  of  moments  may  be  due  to  the 
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fact  that  the  magnetic  field  of  an  electric  current  is  "cross  polarized" 
with  respect  to  the  electric  current.  For  example,  the  magnetic  field 
of  an  electric  current  element  is  always  perpendicular  to  that  current 
element  |7,  F.q.  (2-113)1. 

Although  the  nx  is  vital  to  [14,  Fq .  (1)],  we  decided  not  to 
multiply  (7d)  by  n  '  because  of  the  following  two  reasons. 

1)  The  magnetic  field  in  (7d)  as  it  stands  is  not  all  cross 
polarized.  Only  the  magnetic  field  due  to  the  electric 
currents  J  +1  is  cross  polarized.  The  magnetic  field  due  to 
the  magnetic  current  M  is  not  cross  polarized. 

2)  If  (7d)  were  multiplied  by  n  ,  we  would  have  to  calculate 
matrix  elements  not  only  for  the  operators  F  (.1,0)  and 

II  (.1,0)  but  also  for  n  x  F  (J,0)  and  n  x  M  (J,0).  If  the 
functions  which  multiply  u  and  in  (10)  were  the  same, 
then  the  matrices  for  n  x  F  (.1,0)  and  n  /  H  ( J ,  0 )  could  be 
obtained  by  merely  rearranging  the  elements  of  the  matrices 

i  I 

for  F  (.1,0)  and  II  (J,0).  Unfortunately,  the  function  in 
(1 0b)  is  different  from  that  in  (10a)  so  that  multiplication  of 
(7<1 )  by  n  /  would  indeed  require  calculation  of  additional 
matrix  elements. 


r 

According,  to  (32),  the  solution  of  (31)  for  determines  the 

coefficients!'.,  I  .,  and  V  in  (19)-(21).  Thanks  to  (16),  (17), 
nj  n]  ii| 

and  (22),  expressions  (I9)-(21)  for  the  electric  and  magnetic  currents 
become* 
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Jf  + I  =  I 


M  -1 


,l-,t 


2M  -1 
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(41  ) 
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M-2  2M-3 

v„j4,j«c  +  .  I  , 
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j^+ 

n j— n , j +M+-M+2 


(42) 


IV.  ELIMINATION  OF  T  FOR  NEGATIVE  VALUES  OF  n 
- n - 

At  first  glance,  solution  of  (31)  requires  calculation  of  the 

moment  matrices  T^  for  both  negative  and  non-negative  values  of  n. 

However,  we  will  show  that  each  element  of  T  is  either  even  or  odd 

n  -> 

in  n.  Consequently,  by  decomposing  the  excitation  vector  B  into  its 

n 

even  and  odd  parts  with  respect  to  n,  we  will  be  able  to  express  the 

solutions  to  all  of  equations  (31)  in  terms  of  the  solutions  to  matrix 

equations  involving  only  the  moment  matrices  T  for  the  non-negative 

n 

values  of  n. 


The  equations  and  unknowns  in  (31)  were  ordered  so  that  the  moment 
matrix  could,  as  in  (33),  be  expressed  in  terms  of  a  minimum  number  of 
submatrices.  Unfortunately,  the  arrangement  of  equations  and  unknowns 
in  (31)  is  not  suitable  for  comparing  the  elements  of  T  with  those 

of  Tr.  Hence,  we  rearrange  the  equations  and  unknowns  in  (31)  to  obtain 
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(43) 
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In  (43),  Z  ,  r  being  either  t  or  <}>  and  s  being  either  t  or  (J),  is  the 
n 

submatrix  of  all  the  elements  of 


n  z11' 

r  n 


n  z21- 

r  n 


n  z12' 

r  n 


„22+  ,  22 
Z  +  n  Z 
n  r  n 


for  which  the  testing  function  is  r  type  and  the  expansion  function  is 

s  type.  Expression  (44)  was  drawn  from  the  right-hand  side  of  (33).  The 

testing  functions  W^T  and  expansion  functions  J^7  enter  (44)  by  means 
— ni  — nj 

of  (35).  The  s  type  expansion  functions,  s  being  either  t  or  <}>,  are 

those  on  the  right-hand  sides  of  (16),  (17)  and  (22)  with  the  superscript  s. 

The  r  type  testing  functions,  r  being  either  t  or  <J),  are  the  complex  con- 

r  s 

jugates  of  the  r  type  expansion  functions.  Still  in  (43),  Y^  ,  r  being 
either  t  or  if  and  s  being  either  t  or  <J>,  is  the  submatrix  of  all  the  ele- 


32+  32- 

Y  +  Y  ] 
n  n 


for  which  the  testing  function  is  r  type  and  the  expansion  function  is 

3rs  3  r  s 

s  type.  The  submatrices  Z  and  Y  in  (43)  are  similarly  defined  in 

n  n 


terms  of  the  remaining  submatrices 


[Zn3+  +  <f >  Zn3'l 

r 


-y23+-y23- 

L  n  n  -J 

on  the  right-hand  side  of  (33). 

+3s 

The  column  vector  1^  in  (43),  s  being  either  t  or  <}>,  consists  of 
all  the  elements  of 


(45) 
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for  which  the  superscript  on  the  right-hand  side  of  (18)  is  s  and  all 

the  coefficients  I  .  which  multiply  the  s  type  expansion  functions 

nJ  -3  s 

in  (19).  Expression  (45)  was  drawn  from  (32).  Similarly,  V  in  (43) 

n 

is  the  column  vector  of  all  the  coefficients  of  the  s  type  expansion 
functions  in  (21). 

The  column  vector  V  on  the  right-hand  side  of  (43),  r  being 
either  t  or  *,  consists  of  all  the  elements  of 


;il+ 


-il- 


(46) 


for  which  the  testing  function  is  r  type.  Expression  (46)  was  drawn 

—i3r 

from  (38).  Finallv,  the  column  vector  I  in  (43)  consists  of  all  the 

n 

elements  of  I  for  which  the  testing  function  is  r  type.  The  super¬ 
scripts  t  and  <{)  in  (43)  denote  the  types  of  testing  and  expansion  func¬ 
tions.  However,  the  superscript  3  in  (43)  has  no  special  meaning.  It 
serves  only  to  distinguish  matrices  with  it  from  those  without  it. 

What  happens  to  the  submatrices  Zrs  and  Z^rs  on  the  le't-hand 

n  n 

j*  g 

side  of  (43)  when  n  is  replaced  by  -n?  Previously,  the  eleme:.t':«.  of  Z 

3rs  n 

and  Z  were  traced  back  to  those  of  the  elements  (35)  for  which  the 
n 

testing  functions  are  r  type  and  the  expansion  functions  are  s  type. 
Thanks  to  (16),  (17),  (22),  (23),  (24),  (27),  (28),  (30),  and  (37),  the 
matrix  element  (35)  can  be  written  as 


ZPV  m  _  ^Mri  1_  E\jSt 
nij  ni”  i  '4j” 


(47) 


where  J  are  the  expansion  functions  appearing  in  (13)  and  (14)  and 
“ nj 


r  +  .  r  +  * 

■  <4h> 


(48) 


The  subscripts  and  superscripts  on  the  left-hand  side  of  (47)  indicate 
the  testing  function  W^“  and  the  expansion  function  j .  It  is  under¬ 
stood  from  (35)  that  ZP^"f  indicates  WP*  and  Jqt  and  that  7.pc}~  indicates 

nij  — ni  -nj  nij 

WP  and  .  In  (47),  r  =  t  if  WP‘  is  a  t  type  expansion  function. 
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and  r  =  <t>  if  is  a  <f>  type  expansion  function.  Similarly,  s  in 
(47)  is  determined  by  the  type  of  the  expansion  function  J^~ . 
According  to  (16),  (17),  and  (22),  the  subscript  j'  in  (47)  is 


Similarly,  the  subscript  i'  in  (47)  is 


q  -f 

also  determined  bv  J  . , 

—nj 

determined  by  Explicit  expressions  for  i'  and  j'  are  not 

needed  now.  If  the  testing  functions  and  expansion  functions  in 
[12,  Eqs.  (2)— (5) ]  are  decorated  by  appending  1  to  their  super¬ 
scripts  and  by  changing  their  subscripts  i  and  j  to  i'  and  j'  and 
if  the  normalized  electric  field  operator  -E  /n  used  in  [12]  is 

+  4; 

replaced  by  -E  /n  ,  then  the  elements  of  the  moment  matrix  used  in 
[12]  will  be  coerced  into  becoming  the  right-hand  side  of  (47). 

The  moment  matrix  used  in  [12]  is  given  by  [12,  Eqs.  (9)— (12) ] 
from  which  it  is  evident  that  the  tt  and  (p(p  submatrices  are  even  in  n 
and  that  the  <()t  and  tcj)  submatrices  are  odd  in  n.  Hence,  the  super¬ 
scripted  Z^  submatrices  on  the  left-hand  side  of  (43)  satisfy 
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happens 

to  the 

submatrices 

YrS  and  Y 

(49a) 


(49b) 


3rs 


on  the  left-hand 


side  of  (43)  when  n  is  replaced  by  -n?  Previously,  the  elements  of  Y 

3rs  1 

and  Y  were  traced  back  to  those  of  the  elements  (36)  for  which  the 
n 

testing  functions  are  r  type  and  the  expansion  functions  are  s  type. 

p 1  q  *  r  +  g  * 

Replacement  of  Wr,  and  J  ,  in  (36)  by  W  ”,  and  J  ,,  gives 

-ni  -nj  7  — ni  -nj 


rs 


Ypq:  =  -  v 


nij 


ni 


4-  C  * 

’  » 


0)>, 


(50) 


The  change  in  testing  functions  and  expansion  functions  in  going 
from  (36)  to  (50)  is  the  same  as  the  change  that  sent  (35)  into 
(47). 
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Tf  the  testing  functions  and  expansion  functions  in  [14] 

are  decorated  by  appending  the  superscript  +  and  by  changing  their 

subscripts  i  and  j  to  i'  and  j'  and  if  the  magnetic  field  operator 
s  * 

H  used  in  [14]  is  replaced  by  ii  ,  then  the  elements  of  the  moment 

rs  * 

matrix  in  [14]  denoted  by  Y  become 


where 


ni  J 


<W  n 

— ni  — 


f  c  4 

H  (J  . 

- 


0)> 


-f 


(51) 


(52) 


According  to  (52),  n  is  a  unit  vector  normal  to  the  generating  curve 

of  (S  +A) .  Although  we  have  been  omitting  the  vector  designation  from 

vectors  inside  the  symmetric  product,  we  decided  to  designate  n*  and 
* 

H  as  vectors  in  (51)  to  clearly  indicate  that  the  vector  product 

4  4 

n  ■<  H  is  intended  there.  Striving  toward  (50),  we  rewrite  (51)  as 


r*  c  +  ♦  r*  t  *  <2 

!  =  <n  '  W  H  (J  0)-, 

in  j  -  -ni  — nj 


|r  "  c’* 
)s  =  t  .4* 


(53) 


Modified  testing  functions  W  are  now  defined  bv  inter- 
,  -ni  f , 

changing  the  t  and  components  of 


=  v(C 

•  V 

( 54a) 

=  u.(W*\, 

"Til 

•  V) 

(54b) 

r  *  n  r + 

Replacement  of  W  by  W  ..  in  (53)  gives 

r  -1U  ’Hi' 
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*S'Y  =  *  •  "'<y  *  °> 


js  =  t  ,<p 


The  symbol  ~  on  the  left-hand  side  of  (55)  indicates  that  the  testing 

j*  f  j*  * 

function  is  W  ..  instead  of  W  ...  In  view  of  (52)  and  (54),  (55) 

— ni  -ni 

becomes 


<wni"  H±<-> 


(t'+  ♦  q  t  f-  o  + 

-  <»„!'•  H  <4f’  °>"  ■  Yni'i' 


s  =  t,<p 


The  left-hand  sides  of  (56)  are  the  r  =  t  and  r  =  <P  portions  of  (50). 
Hence,  (56)  states  that  the  r  =  t  matrix  element  in  (50)  is  the  negative 
of  the  r  =  <p  matrix  element  modified  from  [14]  and  that  the  r  =  <j>  matrix 
element  in  (50)  is  the  r  =  t  matrix  element  modified  from  [14].  The 
expression  "matrix  element  modified  from  [14]"  denotes  the  matrix  element 

IT  Aj*  + 

of  [14]  modified  by  replacing  the  testing  function  by  W^.,  of  (54), 

S  s  *■ 

replacing  the  expansion  function  by  J  , ,  and  replacing  the  magnetic 
field  operator  HS  by  Hi  Note  that,  since  all  Ampere's  law  contributions 
were  suppressed  in  (3c)  and  (3d),  the  operator  H'  does  not  contain  any 
Ampere's  law  contribution. 

It  is  evident  from  [14,  Kq.  (12)]  that  the  tt  and  matrix 
elements  on  the  right-hand  sides  of  (56)  are  even  in  n  and  that  the  <J)t 
and  t 0  matrix  elements  on  the  right-hand  sides  of  (56)  are  odd  in  n. 
Consequently,  the  tt  and  '{>4>  matrix  elements  in  (50)  are  odd  in  n  and  the 
ip t  and  t<J)  matrix  elements  in  (50)  are  even  in  n.  Hence,  the  submatrices 

rq  ^rq 

Y  '  and  Y  on  the  left-hand  side  of  (43)  satisfy 
n  n 


(57a) 


(57b) 
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Next,  the  excitation  vector  on  the  right-hand  side  of  (43)  is 
written  as 


The  last  superscript  0  or  (f  on  the  column  vectors  on  the  right-hand  side 
of  (58)  distinguishes  those  column  vectors  from  each  other  and  from  the 
excitation  vector  on  the  left-hand  side  of  (58).  For  a  particular  value  of 
n,  the  solution  to  (43)  is  the  sum  of  the  solutions  obtained  by  replacing 
the  excitation  vector  by  (59a)  and  (59b)  successively.  Hence,  we  write 
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where  Che  column  vectors  on  right-hand  side  of  (60)  satisfy 
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If  n  is  replaced  by  -n  in  (61),  and  if  (49),  (57),  and  (62)  are  substituted 
for  the  -n  submatrices  and  excitation  vectors,  then  comparison  of  the  re¬ 
sulting  matrix  equation  with  (61)  yields 


< 

i™" 

-n 

v3t:" 

-n 

iwi 

-n 

=  [i3t0 

n 

n 

-v3te 

n 

n’,r 

x 

-n 

v3t^ 

-n 

VW1 

-n 

=  (-i3t* 

n 

n 

v3t* 

n 

V  v  ]  (63a) 
n 

-V3^]  (63b) 
n 


The  solution  to  (43)  is  now  given  by  (60)  where  the  column  vectors  for  non¬ 
negative  values  of  n  on  the  right-hand  side  of  (60)  are  obtained  by  solving 
(61)  and  the  column  vectors  for  negative  values  of  n  are  given  by  (63). 


The  result  stated  in  the  previous  sentence  concerning  the  solution 

to  (43)  will  eventually  be  applied  to  the  solution  to  (31).  The  solution 

X  to  (31)  is  written  as 
n 

X  =  X'1  +  Xs  (64) 

n  n  n 

-V 

where,  for  non-negative  values  of  n,  X^  and  X^  are  obtained  by  solving 


-Ml 

T  X' 
n  n 


B 

n 

-Mi  -►.*> 

T  X  =  B’ 
n  u  n 


(65a) 


(65b) 


29 


where 
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The  elements  of  the  column  vectors  on  the  right-hand  sides  of  (67)  are 

P  ^ 

given  by  (39).  In  the  same  way  that  WK",  p-1,2,3,  in  (35)  was  replaced 

— ni 

r  +  n  4 

by  U^,  in  order  to  obtain  (47),  the  testing  function  in  (39)  can  be 
r  + 

replaced  by  W^,.  On  the  right-hand  sides  of  (67),  the  upper  sign  applies 
to  elements  for  which  the  testing  function  in  (39)  is  t  type  and  the  lower 
sign  applies  when  the  testing  function  is  <1  type. 

We  wish  to  express  the  currents  (40)-(42)  in  terms  of  the  elements 

of  X'  and  X  for  non-negative  values  of  n.  With  this  objective  in  mind 
n  n 

and  taking  a  cue  from  (64),  we  separate  (32)  into 
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la  view  of  (68),  the  currents  (40)-(42)  become 
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0  A 
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where  the  currents  with  the  superscript  0  on  the  right-hand  sides  of  (69) — 
(71)  are  due  to  the  set  of  all  the  X^'s  for  all  values  of  n.  Similarly, 
the  currents  with  the  superscript  <p  are  due  to  the  X^’s.  More  explicitly, 

(J*11  +  JU)*  u  being  either  0  or  <t> ,  is  given  hy  (40)  with  I1+  and  I  .  re- 

nj  nj 

placed  by  I  u  and  1^  ,  respectively.  Here,  I1!0  and  IU.  are  the  jth 
nJ  03  nj  nj 

elements  of  1^  and  l“,  respectively.  Likewise,  (J_u  +  Ju)  is  given  by 

(41)  with  I*  and  I  .  replaced  by  the  jth  elements  I^7U  and  IU.  of  I1_u 
nJ  nj  nj  nj  n 

a°d  I  •  Finally,  M  is  given  by  (42)  with  V  .  replaced  by  the  jth  element 

vu.  Of  vu.  nJ 

nj  n 

In  view  of  (63) , 
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where  the  upper  sign  applies  to  the  coefficients  of  the  t  type  expansion 
functions  in  (40)-(42)  and  the  lower  sign  applies  to  the  coefficients  of 
the  <p  type  expansion  functions  in  (40)-(42).  Thanks  to  (72)  and  (10),  the 

i 

t  and  <i>  components  of  the  currents  on  the  right-hand  sides  of  (69)-(7l) 
become 
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where 
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where 


n  =  0 


n  >  0 


In  (73)-(R4),  k  is  the  propagation  constant  in  the  medium  characterized  by 

(li+,  ‘  +).  At  present,  placement  of  the  factor  k+  in  (73)  —  (78)  and  the  factor 

l/k+  in  (79)— (84)  can  only  be  viewed  as  a  matter  of  choice.  The  j  that 

appears  in  the  factor  2j  which  alternates  with  e  in  (79)— (84)  is  /-T.  The 

n 

rest  of  the  j's  in  ( 7 9 )  —  ( 84 )  are  summation  indices. 

Because  of  the  subsectional  nature  of  the  triangle  functions,  only  one 

-f  —-4* 

of  them,  namely  T.  ^(t  ),  contributes  to  (79a)  when  t  =  tj  where  j-1  is  one 
possible  value  of  any  of  the  summation  indices  in  (79a).  Similar  results 

+  -4- 

can  be  stated  for  (80a)  and  (81a).  The  pulse  function  P”(t”)  is  centered 
+  +  J 

at  t  =  t .  where 


vi  ■  \  +  Ei«) 


Because  of  the  subsectional  nature  of  the  pulse  functions,  only  one  of 


them  contributes  to  (79b)  when  t 


tj  where  j  is  one  possible  value  of 


any  of  the  summation  indices  in  (79b).  Similar  results  can  be  stated 
for  (80b)  and  (81b).  Hence,  the  values  of  (79)— (81)  at  the  peaks  of  the 
triangle  functions  and  at  the  centers  of  the  pulse  functions  are  given  by 
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The  mm-7oro  values  of  J 


t+'' 


-+ 


stated  at  t^  and  t^  in  (87a)  are  due  not 

to  peaks  of  triangle  functions  but  to  the  fact  that  tT  and  t*+  lie  on 

1  pp 

the  ?:  axis.  These  non-zero  values  were  established  in  the  paragraph 
which  follows  the  one  containing  (12).  The  starting  and  finishing  values 
of  zero  were  included  in  (89n)  to  emphasize  the  fact  that  goes  to  zero 
at  the  edges  of  the  aperture. 


The  values  of  (82) -(84)  at  the  peaks  of  the  triangle  functions  and 
at  the  centers  of  the  pulse  functions  are  not  written  here.  It  suffices 
to  say  that  these  values  are  given  by  t he  right-hand  sides  of  (87)-(89) 
with  the  (actors  :  and  2j  interchanged  and  with  the  superscript  0  replaced 
by  T. 


In  summary,  the  electric  and  magnetic  currents  are  given  by  (69)-(71) 

and  (73)— (84)  of  which  (79)-(81)  specialize  to  (87)— (89) .  According  to 

(68),  t he  t's  and  V's  in  (79)-(84)  are  elements  of  X  and  .  Now,  X  and 
, .  n  n  u 

Xn  are  obtained  by  solving  (65)  for  non-negative  values  of  n.  The  moment 

matrices  T  in  (65)  are  given  by  (33)— (36) .  However,  (33)  was  reduced  to 

(47)  and  (36)  was  reduced  to  the  combination  of  (50)  and  (56).  The  exci- 

tat  ion  vectors  B  and  V  in  (65)  are  £iven  bv  (66).  (67),  and  (39)  or 
n  4  n 

rather  (  39)  with  W1’ . ,  p  =  1,2,3,  replaced  by  Wr .  , . 
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V.  PLANK  WAVE  L NCI PENCK 

In  this  section,  the  method  of  moments  solution  summarized  in 
the  last  paragraph  of  Section  IV  is  specialized  to  the  case  in  which 
(l’1+,  H1+)  is  an  obliquely  incident  plane  wave  and  R1  =  H1  =0.  If 
the  excitation  (R1+,  H1+)  is  rotated  in  ^  by  an  angle  (J)  ,  then  it  is 
evident  from  rotational  symmetry  that  the  response  will  rotate  by  the 
same  angle  ij  .  Hence,  we  may,  without  loss  of  generality,  assume  that 
the  propagation  vector  of  the  incident  plane  wave  lies  in  the  xz  plane. 
Any  plane  wave  can  be  written  as  the  sum  of  a  0-polarized  plane  wave  and 
a  i-polarized  plane  wave.  Thus,  it  suffices  to  solve  the  electromagnetic 
problem  first  for  a  0-polar ized  incident  plane  wave  and  then  for  a  /p- 
polarized  incident  plane  wave. 

For  the  0-polarized  incident  plane  wave,  we  choose 
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The  1 -polar ized  incident  plane  wave  is  defined  by 
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In  (90)  .and  (91),  r  is  the  radius  vector  from  the  origin,  and  k  is  the 
propagation  vector  given  by 


k  =  -  k  (u  sin 
-t  at 


+  u  cos  n  ) 
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Also, 
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are  unit  vectors  given  by 


(92) 
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In  (92)  and  (93),  u  ,  u  ,  and  u  are  unit  vectors  in  the  x,y,  and  z 

-x  -  -y  — z 

directions,  respectively.  The  "t"  in  (92)  and  (93)  stands  for  transmitter. 
Presumably,  each  of  the  incident  waves  (90)  and  (91)  is  produced  by  a 
distant  transmitter.  According  to  (92),  the  transmitter  bearing  is  given 
by  ( ■•,*,)  =  (■•  ,  0)  where  0  is  the  angle  that  the  radius  vector  from  the 
origin  makes  with  the  positive  z  axis.  When  (6,1))  =  (9  ,  0),  u^  and  u^ 
coincide  with  the  unit  vectors  in  the  0  and  <t>  directions,  respectively. 

The  excitation  vectors  B  and  of  (66)  are  now  obtained  for  the 

n  n  +  r+ 

••-polarized  wave  (90).  Replacement  of  W  ,  p  =  1,2,3,  by  W  in  (39) 

-n i  .  -Til .  .  _ 

and  substitution  of  (90)  and  (0,0)  for  (jy1  ,  H1  )  and  (E1  ,  H1  )  in 
(39)  yield 
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The  values  of  r  and  i1  in  (94a)  depend  on  W  , ,  the  values  of  r  and  i'  in 

2+  ~tl1  3+ 

(94c)  depend  on  W  , ,  and  the  values  of  r  and  i'  in  (94d)  depend  on  W  .. 

T11  M  ij)  +  ’-nl 

It  is  evident  from  (90c)  and  (91b)  that  H  =  -  E  /0  so  that  (94)  becomes 
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Now,  (95)  reduces  to 
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(96a) 
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where  the  V's  on  the  right-hand  sides  of  (96)  are  given  by  [12,  Eqs.  (113), 
(114),  (116),  and  (117)]  with  i  replaced  by  i',  k  replaced  by  k+,  and  the 
quantities  in  the  integrands  made  appropriate  to  the  generating  curve  of 
(S  +A) .  The  reader  should  not  be  misled  into  believing  that  the  right- 
hand  sides  of  (96a)  and  (96c)  are  equal  to  each  other.  Because  r  and  i' 
originated  in  (94),  their  values  in  (96c)  are  not  necessarily  the  same  as 
those  in  (96a). 

With  regard  to  the  right-hand  sides  of  (96),  it  is  evident  from 

[12]  that  Vt0,  and  V44 ,  are  even  in  n  and  that  V4®,  and  Vt4>.  are  odd  in 

ni  ni'  ni  ni' 

n.  Hence,  substitution  of  (96)  into  (67)  with  due  regard  to  (66)  gives 
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B4  =  0  -  (97b) 

n 


Of  course,  the  elements  of  the  row  vectors  on  the  right-hand  side  of 
(97a)  are  given  by  (96). 

The  excitation  vectors  B^  and  B4  of  (66)  are  now  obtained  for  the 

n  n 

([•-polarized  wave  (91).  Replacement  o  f  ^ni»  P  =  i*2’3’  fey  W^,  in  (39), 
substitution  of  (91)  and  (0,  0)  for  (Ei+,  Hi+)  and  (E*  ,  H*  )  in  (39), 
and  use  of  the  relationship  H4  =  ^/n+  yield 
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Now,  (98)  reduces  to 
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where  the  V's  on  the  right-hand  sides  of  (99)  are  given  by  [12,  Eqs.  (113), 

(114),  (11  (3) ,  and  (117)]  with  i  replaced  by  i'  and  k  replaced  by  k+.  Since 

Vt'.  ,  and  V^f,  are  even  in  n  and  V^t  and  V*”?,  are  odd  in  n,  substitution 
ni  ni  ni  ni 

of  (99)  into  (67)  with  due  regard  to  (66)  gives 
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where  the  elements  of  the  row  vectors  on  the  right-hand  side  of  (100b) 
arc  given  by  (99).  The  superscripts  0  and  i}>  that  were  first  introduced 
in  (64)  and  (65)  had  no  mnemonic  meaning  then,  but  seem  to  denote  the  0 
and  ^  polarized  incident  plane  waves  now  in  (97)  and  (100). 


To  avoid  confusion 
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to  achieve  harmony  with  (66) 
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According  to  (99), 
given  by 


: he  it h 
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The  V's  on  the  right-hand  sides  of  (103)  and  (104)  are  given  by  [12,  Eqs. 
(113),  (114),  (116),  and  (117)]  modified  as  indicated  below  (96).  The 
auxiliary  equations  involving  the  W's  in  (103)  and  (104)  determine  r  and 
i'.  For  instance,  the  values  of  r  and  i'  in  (103a)  are  such  that 

wr+,  -  u1*. 

— n  i  — n  l 

It  is  now  evident  that  the  electric  and  magnetic  currents  due  to  the 
F-polarized  incident  plane  wave  (90)  are  given  by  (73)— (75)  and  (79)-(81) 

where  the  I's  and  V's  in  (7 9)  —  (81 )  are  the  elements  (68a)  of  XJ  obtained  by 

-0  n 

solving  (65a)  for  non-negative  values  of  n  with  B  given  by  (101).  Similarly, 

n 

the  currents  due  to  the  fi-polar  ized  wave  (91)  are  given  by  (76)-(78)  and 

(82)-(84)  where  the  I's  and  V's  in  ( 82) — (84 )  are  the  elements  (68b)  of 

.  n 


obtained  by  solving  (65b)  with  given  by  (102).  Because  the  testing 
functions  Wr^,  and  normalized  fields  E^/n+  and  Er/n+  all  have  the  dimension 

~ni  —  - 

of  (1/distance),  the  elements  of  the  excitation  vectors  B  of  (101)  and  IT 

n  n 

of  (102)  are  dimensionless.  Since  the  elements  of  the  moment  matrix  T  are 

0  d)  ^ 

dimensionless,  the  elements  of  X  and  X9  obtained  by  solving  (65)  are  also 

n  n 

dimensionless.  As  a  result,  the  Fourier  coefficients  (79) — (84)  of  the  cur¬ 
rents  are  dimensionless.  This  is  not  surprising  because,  due  to  the  facts 
that  an  electric  current  has  the  same  dimension  as  its  magnetic  field  and 
that  a  magnetic  current  has  the  same  dimension  as  its  electric  field,  the 
electric  currents  (73)— (74)  and  (76)-(77)  must  have  the  same  dimension  as 
(90c)  and  (91c)  and  the  magnetic  currents  (75)  and  (78)  must  have  the  same 
dimension  as  (90b)  and  (91b). 
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VI.  NUMERICAL  RESULTS 

A  computer  program  was  written  to  calculate  the  Fourier  coef¬ 
ficients  (87)-(89)  of  the  electric  and  magnetic  currents  due  to  the 
O-polarized  wave  (90),  which  is  generally  obliquely  incident.  This 
program  will  be  described  and  listed  along  with  sample  input  and  output 
data  in  a  forthcoming  report.  The  'f-polar ized  wave  (91)  was  avoided  for 
fear  of  making  the  program  unwieldy.  According  to  the  exposition  in 
Sections  IV  and  V,  the  Fourier  coefficients  of  the  currents  due  to  the 
^-polarized  wave  can  easily  be  obtained  by  making  two  minor  changes  in 
the  computational  procedure  for  the  9-polarized  wave.  These  changes  are 
replacement  of  (101)  by  (102)  and  interchange  of  the  factors  and  2j 
in  (79)— (81 )  and  (87)-(89). 

In  this  section,  the  magnitudes  of  the  Fourier  coefficients 
(87)-(89)  of  the  electric  and  magnetic  currents  due  to  an  axially  inci¬ 
dent  plane  wave  traveling  in  the  minus  z  direction  are  plotted  for  four 
different  objects  of  revolution.  Specifically,  the  incident  wave  is 
given  by  (90)  with  0  =  0.  At  6^  =  0,  the  excitation  vector  (101),  being 

obtained  from  (12,  Eqs.  (113),  (114),  (116),  and  (117)],  is  zero  for  n  t  ±1. 
Hence,  the  only  non-zero  Fourier  coefficients  in  (87)-(89)  are  those  for 
which  n=l.  These  Fourier  coefficients  were  calculated  by  using  the  com¬ 
puter  program  introduced  in  the  previous  paragraph.  It  is  evident  from 
(73)— (75)  and  (90)  that  expressions  (87)  and  (88)  are  ratios  of  the  elec¬ 
tric  currents  to  the  incident  magnetic  field  and  that  expressions  (89) 
are  ratios  of  the  magnetic  currents  to  the  incident  electric  field. 

In  Figs.  7-9,  the  magnitudes  of  the  Fourier  coefficients  (87)-(89) 
of  the  currents  on  the  object  of  Fig.  6  are  plotted  by  means  of  symbols. 

The  key  for  the  symbols  is  given  by 
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Fig.  6.  Perfectly  conducting  arrow  with  thick  washer  dielectric  region 
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Fig.  8.  Currents  on  object  of  Fig.  6.  Aperture  formulation,  a=0.2\,  e  =1 . 
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Fig.  9.  Currents  on  object  of  Fig.  6.  Aperture  formulation,  a=0.2'-,  j.  *4. 
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*  K"a5| 

in  the  aperture 

(105e) 

O  !«f(tj)l 

in  the  aperture 

( 1 05  f ) 

In  (105a)  and  (105b),  the  plus  sign  applies  on  S+  and  the  minus  sign 
applies  on  S  .  The  surfaces  S+  and  S  merge  on  the  infinitely  thin  con¬ 
ducting  flange  from  C  to  H  in  Fig.  6.  There,  (105a)  and  (105b)  are 
replaced  by 

|j^+°(t^)  -  (t?)|  on  the  thin  flange  (106a) 

|jf+'(tT)  -  ^(t  )|  on  the  thin  flange  (106b) 

1  j  J-  m 

where  «,  is  such  that  both  t+  and  t p  are  located  at  the  same  point  on  the 

1  + 

flange  and  m  is  such  that  both  t.  and  t  are  located  at  the  same  point  on  the 

J  111 

flange.  The  expressions  whose  magnitudes  are  being  taken  in  (106)  are  the 
Fourier  coefficients  of  the  electric  current  induced  on  the  flange.  In 
general,  the  electric  current  induced  on  an  infinitely  thin  perfectly  con¬ 
ducting  surface  is  the  sum  of  the  exterior  and  interior  currents.  However, 
because  -J  is  used  in  Fig.  3,  the  electric  current  induced  on  the  flange 
in  Fig.  6  is  J  -J  .  Thus,  the  difference  rather  than  the  sum  of  coef¬ 
ficients  is  required  in  (106).  All  of  expressions  (105)  and  (106)  are 
plotted  in  each  of  Figs.  8  and  9.  However,  only  the  expressions  on  the 
conducting  surfaces,  namely  (105a),  (105b),  (106a),  and  (106b)  are 
plotted  in  Fig.  7.  Actually,  the  exrressions  plotted  in  Fig.  7  are  not 
(105a),  (105b),  (106a),  and  (106b)  themselves  but  expressions  comparable 
to  them. 

For  simplicity,  the  vertical  axes  in  Figs.  7-9  are  labeled  [j  j . 

+  ^ 

The  horizontal  axes  are  labeled  t  rather  than  t  or  t  .  Here,  t  represents 
arc  length  between  the  points  A,B,C,...  designated  in  Fig.  6.  The  t  com¬ 
ponents  of  (the  Fourier  coefficients  of)  the  currents,  namely  (105a),  (105c), 
(105e),  and  (106a),  are  plotted  at  the  tick  marks  in  Fig.  6.  These  tick 
marks  represent  the  data  points  (11).  The  $  components  of  the  currents, 
namely  (105b),  (105d),  (105f),  and  (106b),  are  plotted  at  points  midway 
between  the  tick  marks  in  Fig.  6.  In  Figs.  7-9,  the  symbols  used  to 
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plot  (105)  and  (106)  are  connected  by  straight  line  segments  to  improve 
readability.  To  avoid  congestion,  the  aperture  region  from  H  to  F  in 
Fig.  6  is  depicted  twice  in  each  of  Figs.  8  and  9,  once  for  plotting  the 
electric  currents  (105c)  and  (105d)  and  once  for  plotting  the  magnetic 
currents  (105e)  and  ( 1 05 f ) .  The  following  argument  reveals  that  the  t 
axes  in  Figs.  7-9  are  not  always  drawn  to  scale.  Not  all  data  points 
in  Fig.  6  are  equally  spaced.  For  instance,  the  data  points  from  F  to 
G  are  not  quite  as  dense  as  those  from  A  to  B.  More  obviously,  the 
data  points  in  the  aperture  are  twice  as  dense  as  those  from  A  to  B. 
However,  all  data  points  are  equally  spaced  along  the  t  axes  in  Figs.  7-9. 
Incidentally,  the  proper  density  of  data  points,  or  rather  triangle 
functions,  in  the  aperture  is  discussed  in  [8]. 


The  currents  in  Fig.  7  were  obtained  by  means  of  the  E-field  solu¬ 
tion,  the  currents  in  Figs.  8  and  9  by  means  of  the  aperture  formulation, 
and,  if  the  reader wi  1 1  please  glimpse  forward,  the  currents  in  Fig.  15  by 
means  of  the  H-field  solution.  The  E-field  solution  is  that  of  [12]  with 


nt  =  2 

n.  =  20 

9 

nT  =  2 


(107) 


where  n^,  n^,  and  n^  appear  in  the  Gaussian  quadrature  formulas  [12, 

Eqs.  (62),  (64),  and  (132)],  respectively.  The  H-field  solution  is  that 
of  [14]  with  the  values  of  nt,  n^,  and  therein  given  by  (107).  Speci¬ 
fically,  nfc,  n^,  and  nT  appear  in  [14,  Eqs.  (35),  (36),  and  (79)].  The 
roles  played  by  n^,  n^ ,  and  n^  in  [14]  are  similar  to  the  roles  played  by 
nt>  n^,  and  n^  in  [12],  The  aperture  formulation  refers  to  the  solution 
(87)-(89)  developed  in  the  present  report.  In  view  of  (68a),  the  I's  and 
V's  in  (87 ) — (89)  are  given  by  the  solution  X®  to  (65a).  In  Section  IV, 
the  elements  of  the  moment  matrix  T^  were  expressed  in  terms  of  the  ele¬ 
ments  of  matrices  similar  to  the  moment  matrices  used  in  [12]  and  [14]. 

■*■0 

In  Section  V,  the  elements  of  in  (65a)  were  expressed  in  terms  of  the 
elements  of  the  plane  wave  excitation  vectors  used  in  [12].  The  currents 
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attributed  to  the  aperture  formulation  were  obtained  by  calculating 
and  with  nf ,  n^,  and  given  by  (107). 

The  E-field  solution  for  the  electric  current  induced  on  the 
conducting  surfaces  of  the  object  in  Fig.  6  was  constructed  by  using 
two  generating  curves.  The  first  one  is  ABCDEFG  and  the  second  one 
extends  from  the  tick  mark  directly  below  C  to  the  point  H.  Since  the 
computer  program  in  [12]  is  designed  to  handle  only  one  generating  curve, 
the  two  (short)  curves  were  connected  to  form  one  long  curve,  and  then  all 
matrix  elements  associated  with  the  line  connecting  the  two  short  curves 
were  deleted.  The  matrix  elements  associated  withthefirst  pulse  func¬ 
tion  on  the  second  short  curve  were  also  deleted  because  that  pulse  func¬ 
tion  is  a  duplicate  of  one  of  the  pulse  functions  on  the  first  short  curve. 
With  this  arrangement,  two  independent  triangle  functions  are  centered  at 
C.  Because  each  of  these  triangle  functions  satisfies  Kirchhoff's  current 
law  at  C,  the  resulting  E-field  solution  for  the  t  component  of  the  induced 
electric  current  satisfies  Kirchhoff's  current  law  at  C.  The  t  component 
of  the  induced  electric  current  obtained  by  means  of  the  aperture  formu¬ 
lation  also  satisfies  Kirchhof f ' s ' fcurrent  law  at  C.  However,  satisfaction 
of  Kirchhoff's  current  law  by  the  t  component  of  the  induced  electric 
current  at  C  is  masked  in  Figs.  7-9  because  magnitudes  of  currents  rather 
than  real  and  imaginary  parts  of  currents  are  plotted  there. 

■^he  t  component  of  the  induced  electric  current  is  zero  at  H  in 
Fig.  7  because  no  triangle  function  is  centered  there.  The  t  component 
(106a)  of  the  induced  electric  current  is  zero  at  H  in  Figs.  8  and  9  be¬ 
cause  the  two  J's  in  (106a)  are  equal  to  each  other  at  H.  These  two  J's 
have  to  be  equal  to  each  other  at  H  because  they  are  both  equal  to  the 
aperture  current  there. 

It  is  obvious  that  the  E-field  solution  for  the  t  component  of 
the  induced  electric  current  has  to  be  continuous  at  F  in  Fig.  6.  The 
following  reasoning  shows  that  the  t  component  of  the  induced  electric 
current  obtained  by  using  the  aperture  formulation  is  also  continuous  at  F. 
Because  the  t  direction  is  the  same  as  the  t+  direction  on  the  G  side  of  F, 
the t component  of  the  induced  electric  current  is  (t  )  there.  Now, 
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t  and  t  are  oppositely  directed  on  the  E  side  of  F  and  there  is  a 
compensating  minus  sign  on  in  Fig.  3  so  that  the  t  component  of  the 
induced  electric  current  is  .l|  ^(t  )  on  the  E  side  of  F.  However, 

,J^+^(t+)  =  ^ ( t  )  at  F  because  they  are  both  equal  to  the  aperture 

current  there.  Hence,  the  t  component  of  the  induced  electric  current 
obtained  by  using  the  aperture  formulation  is  continuous  at  F. 

The  notation  previously  developed  with  specific  reference  to 
Figs.  6-9  and  15  applies  not  only  to  these  figures  but  also  to  Figs. 

10-14  and  16-22.  The  object  in  Fig.  10  differs  from  the  exemplary 
object  in  Fig.  1  in  that  neither  the  exterior  surface  S+  nor  the  in¬ 
terior  surface  S  has  a  part  above  the  aperture  in  Fig.  10.  Nevertheless, 
the  computer  program  mentioned  in  the  first  paragraph  of  this  section  is 
still  applicable.  Figure  11  shows  currents  on  the  spherical  shell  and 
associated  aperture  region  of  Fig.  10  when  C  =  1.  The  currents  plotted 
from  A  to  B  in  Fig.  11  are  the  electric  currents  (106)  induced  on  the 
conducting  shell. 

The  following  reasoning  shows  that  the  E-field  solution  for  the 
electric  current  induced  on  the  conducting  shell  in  Fig.  10  is  exactly 
the  same  as  the  solution  obtained  by  means  of  the  aperture  formulation 
when  e  =  1.  When  partitioned,  (31)  consists  of  four  matrix  equations. 
The  first  two  of  these  equations  are 


z11+  f1+  +  z12+  I  -  y13+  V  =  vil+ 

n  n  n  n  n  n  n 


(108a) 


11-  -»1-  IT-  v  13-  -v  ‘■il- 

:|Z  I  +  n  7.  I  -  V  =  v“  (108b) 

rnnrnnnnn 


Since  the  conducting  surface  in  Fig.  10  is  infinitely  thin,  the  surfaces 
S+  and  S  exemplified  in  Fig.  1  are  identical.  Furthermore,  f  =  1  and 
0  =  1  so  that  the  superscripts  +  and  -  attached  to  the  7,'s  and  Y's  in 
(108)  lose  their  distinction.  Hence,  the  difference  between  (108a)  and 
(108b)  reduces  to 


Z11+[T1+  -  I1"]  =  Vil+ 
n  n  n  n 


ti  1  - 


(109) 


Because  the  elements  of  [1 


1  + 


rl-. 


can  be  interpreted  as  the  unknown 


Fig.  10.  Infinitely  thin  perfectly  conducting  spherical  shell  with 


Currents  on  object  of  Fig.  10.  Aperture  formulation, 
k+a  «  2.5,  c  -  1. 


Fig.  11. 
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coefficients  appearing  in  the  expansion  of  the  electric  current  induced 
on  the  shell,  (109)  is  precisely  the  matrix  equation  for  the  E-field  solu¬ 
tion.  As  a  result,  the  E-field  solution  for  the  electric  current  in¬ 
duced  on  the  shell  in  Fig.  10  is  exactly  the  same  as  the  electric  current 
plotted  from  A  to  B  in  Fig.  11. 

The  magnetic  currents  (105e)  and  (105f)  plotted  from  B  to  C  in 
Fig.  11  are  comparable,  respectively,  to  the  functions  jg(h)|  and  |f(0)j 
appearing  in  [15,  Fig.  7].  The  horizontal  axes  in  [15,  Fig.  7]  and  Fig.  11 
run  in  opposite  directions.  The  statement  that  ka  =  2.75  in  the  caption 
of  [15,  Fig.  7]  is  incorrect  because  ka  is  really  2.5  there.  The  currents 
shown  in  Fig.  12a  were  calculated  by  using  the  data  points  designated  in 
Fig.  10.  However,  the  currents  shown  in  Fig.  12b  were  calculated  with 
the  curve  from  A  to  B  in  Fig.  10  divided  into  not  10  but  15  equal  intervals. 
For  Fig.  12b,  the  data  points  in  the  aperture  were  those  designated  in 
Fig.  10. 

The  object  in  Fig.  13  was  obtained  by  taking  a  perfectly  conducting 

solid  sphere  of  radius  a  centered  at  the  origin  and  replacing  the  conducting 

material  in  the  region  |z|  <  a  sin  15°  by  dielectric  material  characterized 

bv  (E  ,  e  e  ).  The  E-field  and  H-field  solutions  for  the  electric  current 
o  r  o 

induced  on  the  conducting  surfaces  of  Fig.  13  are  shown  in  Figs.  14  and  15, 
respectively.  These  solutions  were  obtained  by  connecting  the  generating 
curves  ABC  and  DEF  in  Fig.  13  to  obtain  the  single  curve  ABCDEF  and  then 
deleting  the  matrix  elements  associated  with  the  connecting  line  CD.  Di¬ 
visions  by  zero  were  avoided  by  arbitrarily  setting  k+p  equal  to  1  midway 
between  C  and  D.  Here,  p  is  the  distance  from  the  z  axis. 

The  currents  shown  in  Figs.  16  and  17  are  on  the  object  of  Fig. 

13  and  were  obtained  from  the  aperture  formulation.  The  object  in  Fig. 

13  differs  from  the  exemplary  object  in  Fig.  1  in  that  the  parts  of  S 
above  and  below  the  aperture  are  disjoint  in  Fig.  13.  However,  the  com¬ 
puter  program  mentioned  in  the  first  paragraph  of  this  section  allows  for 
this.  The  induced  electric  currents  from  A  to  C  and  from  D  to  F  are 
slightly  different  from  each  other  in  Figs.  14-16.  The  induced  electric 
current  on  the  conducting  surfaces  of  Fig.  13  obtained  from  the  aperture 

formulation  when  c  “1  differs  from  the  E-field  solution  for  the  induced 
r 


Currents  on  object  of  Fig.  13.  Aperture  formulation, 
k+a  =  2.5,  »  2. 
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electric  current  because  of  the  following  reasoning.  Of  the  triangle 
functions  used  in  the  aperture  formulation,  those  centered  at  B  and  E 
are  essential  to  the  induced  electric  current.  However,  these  triangle 
functions  extend  into  the  aperture.  On  the  other  hand,  all  of  the  triangle 
functions  used  in  the  E-field  solution  are  confined  to  the  conducting  sur¬ 
faces.  In  this  respect,  the  aperture  formulation  differs  from  the  E-field 
solut ion. 

Figures  19-22  show  currents  on  the  object  of  Fig.  18.  On  the  object 
of  Fig.  18,  EFB  is  the  generating  curve  of  the  aperture.  The  object  in 
Fig.  18  differs  from  the  exemplary  object  in  Fig.  1  in  that  neither  the  ex¬ 
ternal  surface  S+  nor  the  internal  surface  S  has  a  part  below  the  aperture 
in  Fig.  18.  However,  the  computer  program  mentioned  in  the  first  paragraph 
of  this  section  allows  for  this.  As  seen  from  Figs.  19-21,  the  E-field 
solution,  the  H-field  solution,  and  the  aperture  formulation  each  give 
slightly  different  results  for  the  electric  current  induced  on  ABCD  in 
Fig.  18.  The  E-field  solution  and  the  aperture  formulation  give  different 
results  for  this  current  because  in  the  aperture  formulation  the  conductor 
is  coupled  to  the  aperture  at  point  B  in  Fig.  18. 
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